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A. LISTINGS OF MAXLABEL PROGRAMS A-1 
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1 . SCOPE 


This document describes the program to represent the CLASSY cluster 
distributions as equi probability ellipses In brightness/greenness space. 

The program Is written exclusively In Fortran IV-G. 
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2. APPLICABLE DOCUMENTS 


The following document forms a part of this specification: 
TIRF 80-0012 Graphical Output of CLASSY Distribution Clusters. 


2-1 




3. SYSTEM DESCRIPTION 


The Equlprobability Ellipses Program was designed and implemented to provide 
graphics output so that the cluster distributions estimated by CLASSY may 
be visualized. The approach is to project the hyper-ellipsoid for each 
cluster distribution onto the brightness greenness plane for each acquisi- 
tion. This projection is done in two ways. The first projection 
corresponds to the selected brightness-greenness channels with the other 
channels unspecified. The second projection corresponds to a projection 
onto the selected brightness greenness channels with the other channels 
fixed at the mean value. Each of these projections produce an ellipse with 
the ellipse for the second projection contained in the ellipse for the first 
projection. 


The equations for the two ellipses plotted for each cluster are as 
follows; 


(1) (X, - ;! * 2(x, 


v.,)(Xj 


“j' “ )! * '”1 ■ if 


XX D 


(2) (X. - + 2(X. - p.)(X^^ - p.) 0 ^^ + {X. - = D 


Where 


-1 , m 
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0. .0 . . a. . 
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n JJ 


o . . 
- TJ 


'^ij ~ element of the cluster covariance matrix. 

Note that o’’ and o’"^ are elements of the matrix which is the inverse 
of the cluster covariance matrix. 
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D Is a value to be selected by the user. The values for the mean vector, 
and Inverse covariance matrix should be read from a data file produced by 
CLASSY. 

The user may elect to plot one or both of the ellipses In the ellipse pair. 
The ellipse represented by equations (1) and (2) may be plotted two ways; 

(1) The program will plot the ellipse or ellipse pair for each cluster 
for a given acquisition on a single page using the same symbol for 
each ellipse pair for a given cluster. Each plot has the same scale. 

(2) The program will plot the ellipse or ellipse pair for each acquisi- 
tion for a given cluster on four pages using the same symbol for each 
ellipse pair for a given acquisition. The user must reconstruct the 
4 pages of plot data as one 2 page by 2 page plot. Each plot has 

the same scale. 

The capability Is Implemented as a stand alone program with enough modularity 
so that It might be easily modified In the future. The output Is 
generated on the line printer with user provided headings and axes for each 
plot. The program has the capability of handling up to 30 clusters 
and 6 acquisitions. 

The driver program for the clustering system Is Ellipses. The calculation 
subprograms are ACQ1CL to calculate the ellipses for all of the acquisitions 
for one CLASSY cluster and CLIACO to calculate the ellipses for all of the 
CLASSY clusters for one acquisition. 

The plotting subroutines are PLOTPT which stores the points In the plot 
matrix, PLOTCL which prints the plot for one CLASSY cluster and PLOTAC 
which prints the plot for one acquisition. 

CLASSY cluster data Is read from the file by READCC and data for the current 
CLASSY cluster Is stored by 6ETCC. 
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3.1 HARDWARE DESCRIPTION 


The Equiprobablllty Ellipses Program Is operational on the IBM 3031 at 
Purdue LARS under the CMS370 operating system. The program utilizes the 
IBM Fortran IV-G computer. 

3.2 ELLIPSE EXEC FILE 

ELLIPSE Is called by the following EXEC file: 

ELLIPSE (Segment number) 
where segment number > number of se^nt. 


3.3 SOFTWARE DESCRIPTION 


3.3.1 SOFTWARE COMPONENT NO. 1 (ELLIPSES) 

ELLIPSES represents the CLASSY cluster distributions as equi probability 
ellipses In brightness/greenness space. 

3. 3. 1.1 LI nkaqes 

ELLIPSES Is the driver program. ELLIPSES calls ACQ1CL. CL1ACQ and I^ADCC. 

3. 3. 1.2 Interface 

The blank coninon block and /CC/ coninon block are used as Interfaces with the 
subprograms . 

3. 3. 1.3 I:iputs 


TITLE 

Blank Common 

Plot title 

XLABEI. 

Blank Connon 

Label for X axis 

YLABEL 

Blank Connon 

Label for Y axis 

D 

Blank Common 

Computational constant 

STARTX 

Blank Common 

Starting X value 

I GROUP 

Blank Common 

Option for group of ellipses 

I OPT 

Blank Coninon 

Option for outer or Inner ellipses 

MQ 

/CC/ 

Number of channels 

lOUT 

/CC/ 

Printer unit 

ICLUNT 

/CC/ 

CLASSY data unit 


3. 3. 1.4 Output 
Same as Input. 

3. 3. 1.5 Storage Requirements 
Not applicable. 
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3. 3. 1.6 Description 

This program establishes the user requirements for the representation of 
the CLASSY cluster distributions and calls tlw subroutines to calculate the 
points of the ellipse and create plots on the line printer. 

3. 3. 1.7 Flowchart 
Not applicable. 

3. 3. 1.8 Listing 

See Appendix A for program. 
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3.3.2 SOFTWARE COMPONENT NO. 2 (ACQICL) 

ACQ1CL calculates the y values for 60 x values for the acquisitions (sets of 
greenness and brightness values) for one CLASSY cluster and calls subprograms 
to store data in the plot and print the plot. 

3. 3. 2.1 Linkages 

ACQICL is called by ELLIPSES. ACQICL calls GETCC, MINV, PLOTPT and PLOTCL. 


3. 3. 2. 2 Interfaces 


The blank common block and /CC/ common block are used as interface between 
ACQICL and ELLIPSES and the subprograms. A calling sequence is used for 
interface with MINV. 


3.3. 2. 3 

Inputs 


MQ 

/CC/ 

Number of channels 

NOCC 

/CC/ 

Number of CLASSY clusters 

eVkIN 

/CC/ 

Inverse of CLASSY covariance matrix 

CCOVAR 

Calling Sequence 

CLASSY covariance matrix 

CMEANS 

/CC/ 

Means for CLASSY clusters 

D 

Blank Common 

Computational constant 


3. 3. 2. 4 Outputs 

Y1 Blank Common 
Y2 Blank Common 
Z1 Blank Common 
Z2 Blank Common 
LINE Blank Common 


Y value for outer ellipse 

Y value for outer ellipse 

Y value for inner ellipse 

Y value for inner ellipse 
Matrix filled with blanks 


3. 3. 2. 5 Storage Requirements 
Not applicable. 
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3. 3. 2. 6 Descriptio n 

For each CLASSY cluster the following events occur. Each set of greenness 
and brightness values are used to calculate ellipses using the covariance 
matrices and means for the CLASSY cluster using the formulas given in 3,0. 
The subprogram PLOTPT is called to store the points in the plot array, and 
PLOTAC is called to print the plots. 

3. 3. 2. 7 Flowchart 
Not applicable. 

3. 3. 2. 8 Listing 

See Appendix A for program. 
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3.3.3 SOFTWARE COMPONENT NO. 3 (CLIACQ) 

CLIACQ calculates the y values for 60 x values for the CLASSY clusters for one 
acquisition (set of greenness and brightness values) and calls subprograms 
to store data in plot and print plot. 

3. 3. 3.1 Linkages 

CLIACQ Is called by ELLIPSES. CLIACQ calls GETCC, MINV, PLOTPT and PLOTCL. 

3. 3. 3. 2 Interfaces 

The blank common block and /CC/ common block are used as Interface between 
CLIACQ and ELLIPSES and the subprograms. A calling sequence Is used for 
Interface with MINV. 


3. 3. 3. 3 

Inputs 


m 

/CC/ 

Number of channels 

NOCC 

/CC/ 

Nun*er of CLASSY clusters 

CVRIN 

/CC/ 

Inverse of CLASSY covariance matrix 

CCOVAR 

Blank Conmon 

CLASSY covariance matrix 

CMEANS 

/CC/ 

Means for CLASSY clusters 

D 

Blank Common 

Computational constant 

3. 3. 3. 4 

Outputs 


Y1 

Blank Coimion 

Y value for outer ellipse 

Y2 

Blank Common 

Y value for outer ellipse 

Z1 

Blank Common 

Y value for Inner ellipse 

Z2 

Blank Common 

Y value for Inner ellipse 

LINE 

Blank Common 

Matrix filled with blanks 


3. 3. 3. 5 Storage Requirements 
Not applicable 
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3. 3. 3. 6 Description 

For each acquisition the following events occur. Each set of greenness and 
brightness values are used to calculate ellipses using the covariance 
matrices and means for the CLASSY cluster using the formulas given in 3.0. 
The subprogram PLOTPT is called to store the points in the plot array, and 
PLOTAC is called to print the plots. 

3. 3. 3. 7 Flowchart 
Not applicable. 

3. 3. 3. 8 Listing 

See Appendix A for program. 
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3.3.4 SOFTWARE COMPONENT NO. 4 (PLOTPT) 

PLOTPT stores the symbols for the ellipses in the matrix LINE. 


3. 3. 4.1 Linkages 

PLOTPT is called by ACQICL and CLIACQ. 

3. 3. 4. 2 Interface 

The blank common block is used as interface with the calling subprograms. 


3. 3. 4. 3 

Inputs 


I GROUP 

Blank 

Common 

Y1 

Blank 

Common 

Y2 

Blank 

Common 

Z1 

Blank 

Common 

Z2 

Blank 

Common 

LINE 

Blank 

Common 


Plot ellipse indicator 

Y value for outer ellipse 

Y value for outer ellipse 
Z value for inner ellipse 
Z value for inner ellipse 
Matrix for line symbols 


3. 3. 4. 4 Outputs 

LINE Matrix containing symbols for points on plots. 


3. 3. 4. 5 Storage Requirements 
Not applicable. 


3. 3. 4. 6 Description 

This subroutine stores syni)ols for the ellipses in the matrix LINE depending 
on which ellipses are requested. 


3.3.5 SOFTWARE COMPONENT NO. 5 (PLOTAC) 

PLOTAC prints a plot for displaying the CLASSY cluster ellipses for one 
acquisition. 

3. 3. 5.1 Linkages 

PLOTAC Is called by CLIACQ. 

3. 3. 5. 2 Interface 

Interface Is accomplished thrcjgh blank common and /CC/ common blocks. 


3.3.5 .3 

Input 


STARTX 

Blank Common 

Starting value for X axis 

LINE 

Blank Common 

Matrix of plot values 

TITLE 

Blank Common 

Title for plot 

XTITLE 

Blank Common 

X label for plot 

YTITLE 

Blank Common 

Y label for plot 

NOCC 

/CC/ 

Number of CLASSY clusters 

NOCHAN 

/CC/ 

Number of current channel 


3. 3. 5. 4 Output 

Plot written on line printer. 

3. 3. 5. 5 Storage Requirement 
Not applicable. 

3. 3. 5. 6 Description 

PLOTAC prints the plot title, prints the lines of plot data Including the 
Y axis label and prints the label for the X axis for 2 plot halves. 


n 
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3. 3. 5. 7 Flowchart 


Not applicable. 

3. 3. 5. 8 Listing 

See Appendix A for program. 
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3.3.6 SOFTWARE CO^f>ONENT NO. 6 (PLOTCL) 

PLOTCL prints a plot for displaying the acquisition ellipses for one 
CLASSY cluster. 

3. 3. 6.1 Linkages 

PLOTCL Is called by ACQICL. 

3. 3. 6. 2 Interface 

Interface Is accomplished through blank common and /CC/ common blocks. 


3. 3.6.3 

Input 



STARTX 

Blank 

Commen 

Starting value for X axis 

LINE 

Blank 

Common 

Matrix of plot values 

TITLE 

Blank 

Common 

Title for plot 

XTITLE 

Blank 

Common 

X label for plot 

YTITLE 

Blank 

Common 

Y label for plot 

NOCC 

/CC/ 


Number of CLASSY clusters 

NOCHAN 

/CC/ 


Number of current channel 

3. 3.6. 4 

Output 




Plot written on line printer. 

3. 3. 6. 5 Storage Require ’^ent 
Not applicable. 

3. 3. 6. 6 Description 

PLOTCL prints the plot title, prints the lines of plot data Including the Y 
axis label and prints the label for the X axis for 1 plot. 
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3. 3. 6. 7 Flowchart 


Not applicable. 

3. 3. 6. 8 Listing 

See Appendix A for program. 


3.3.7 SOFTWARE COMPONENT NO. 7 (READCC) 
fSADCC reads the CLASSY cluster data. 

3. 3. 7.1 Li nkage 

READCC is called from ELLIPSES. 

3. 3. 7. 2 Interface 

Interface Is acconpllshed through the /CC/ common block. 


3. 3. 7. 3 

Input 


I OUT 

/CC/ 

Printer unit 

ICLUNT 

/CC/ 

Unit for CLASSY data file 

r*) 

/CC/ 

Number of channels 

3. 3. 7.4 

Output 


RLPRP 

/CC/ 

Relative proportion array 

CCVLRT 

/CC/ 

Normalization factor element 

CDCON 

/CC/ 

Normalization factor element 

CCMEANS 

/CC/ 

Cluster means array 

CCVRIN 

/CC/ 

Covariance matrix 

NOCC 

/CC/ 

Number of CLASSY clusters 


3. 3. 7. 5 Storage Requirements 
Not applicable. 

3. 3. 7. 6 Description 

READCC reads the CLASSY cluster values into cownon block /CC/. The nun*er 
of sets of CLASSY cluster information is saved in NOCC. 
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3. 3. 7. 7 Flowchart 
Not applicable. 


3. 3. 7. 8 Listing 

See Appendix A for program. 


3.3.8 SOFTWARE COMPONENT NO. 8 (GETCC) 

Move one set of CLASSY data from arrays to corresponding scalars. 


3. 3.8.1 Linkage 

GETCC Is called from GETCC. 

3. 3. 8. 2 Interface 

Interface Is accomplished through blank common and /CLASSY/ blocks. 


3. 3. 8. 3 Input 

ICC /blank common/ current class 


3. 3. 8. 4 Output 

I^LPRP /blank common/ 
CVOLRT /blank common/ 
CMEANS /blank conmon/ 
CVRIN /blank common/ 
DCON /blank common/ 


Illative proportion 
Normalization factor element 
Classy cluster mean array 
Covariance matrix 
Normalization factor element 


3. 3. 8. 5 Storage Requirement 
Not applicable. 


3. 3. 8. 6 Description 

GETCC moves data for this specified class from the arrays of CLASSY data. 


3. 3. 8. 7 Flowchart 
Not applicable. 


3. 3. 8. 8 Listing 

See Appendix A for program. 


3.3.9 SOnWARE COMPONENT NO. 9 (KINV) 

MINV calcualtes the Inverse and determinant of the Input matrix. 


3.3.9. 1 Linkages 

MINV Is called by ACQ1CL and CL1ACQ. 

3. 3. 9. 2 Interfaces 

Interface Is acccmpllshes through calling arguments and connon block /CC/. 

3. 3. 9. 3 Input 

MQ /CC/ 

B Calling Argument 

C Calling Argument 

3. 3. 9. 4 Output 

A Calling Argument Inverted matrix 

VOL Calling Argument Determinant of C 

3. 3. 9. 5 Storage Requirement 
Not applicable. 

3. 3. 9. 6 Description 

MINV calculates A > Inverse of C and V(M. the determinate of C. 

3. 3. 9. 7 Flowchart 
Not applicable. 

3. 3. 9. 8 Listing 
See Appendix A for program. 


Nunter of channels 
Scratch area 
Matrix to be Inverted 


'10 
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4. OPERATION 


ELLIPSES Is operational on the IBH 3031 computer at LARS, Nest Lafayette 
Indiana. 

The program and EXEC files can be loaded from tape 3956, file 6. 

ELLIPSES Is executed by entering the following comnarKl after signing on the 
coa^iuter system. 

ELLIPSES (Segment Number) 

where Segment Nimber « number of the segment. 

Control Information Is defined by Interactive dialogue with the user. 

Text output Is on the terminal and line printer which are assigned In the 
EXEC. 
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APPENDIX A 


LISTINGS OF MAXLABEL PROGRAMS 
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